Optimization of congested traffic flow in systems with a localized periodic 

inhomogeneity 
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We study traffic flow on roads with a localized periodic inhomogeneity such as traffic signals, 
using a stochastic car-following model. We find that in cases of congestion, traffic flow can be 
optimized by controlling the inhomogeneity's frequency. By studying the wavelength dependence 
of the flux in stop-and-go traffic states, and exploring their stability, we are able to explain the 
optimization process. A general conclusion drawn from this study is, that the fundamental diagram 
of traffic (density flux relation) has to be generalized to include the influence of wavelength on the 
flux, for the stop-and-go traffic. Projecting the generalized fundamental diagram on the density-flux 
plane yields a 2D region, qualitatively similar to that found empirically [B. S. Kerner, Phys. Rev. 
Lett. 81, 3797 (1998)] in synchronized flow. 
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The theory of traffic flow has been a subject of compre- 
hensive study for more than half a century [p]-p8[ due to 
its theoretical and practical importance. Much attention 
has been devoted to characterizing the different states of 
congested traffic — |l0|] , including synchronized flow and 
stop-and-go traffic. 

Optimization of systems with localized inhomo- 
geneities such as traffic signals or entrance ramps was also 
extensively studied in the last decades |l8|-^4|. Conven- 
tional traffic signals and ramp control theories are able 
to optimize such quantities as the total delay time of 
all drivers in the system, assuming a well defined flux- 
density relation known as the 'fundamental diagram' of 
traffic flow. 

This basic concept of a fundamental diagram was in- 
grained in traffic flow theory for many decades. It was 
believed that the density-flux relation can be displayed 
as a single curve or as a combination of two isolated 
curves (see e^. p8|,[7|). Recently, an experimental study 
of Kerner pUq] shows that such a fundamental diagram 
does not exist. Instead, synchronized traffic displays a 
two-dimensional region in the density-flux plane. That 
is, for a given value of density, there exist a range of pos- 
sible flux values. Mechanisms to optimize traffic flow by 
approaching the highest values of this range have not yet 
been suggested. 

In this paper we investigate whether this new insight 
on the nature of traffic flow can be applied for optimizing 
the flux close to its maximal value for a given congested 
density. For this purpose, systems with periodic localized 
inhomogeneities are studied, using a recent car-following 
model with multiple stable and metastable states p6| . 
Two types of periodic inhomogeneities are considered: 



(a) A signalized intersection, and 

(b) an on-ramp with a signalized entrance. 

Our study focus on cases of oversaturation, i.e. on cases 
where traffic is congested upstream to these inhomo- 
geneities. 

In the common traffic flow microscopic models such as 
in Refs. |lo| [li]], the acceleration a of a car depends on 
its headway Ax, velocity v, and velocity difference Av 
with the car ahead, i.e. a — a(Ax,v, Av). In particular, 
for the model in 
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where T is the safety time gap, D and Ax{) = vT + D are 
the minimal and the optimal distance to the car ahead. 
A, k, and constants, and the function Z is de- 

fined as Z(x) = (x + \x\)/2. In the following numerical 
solutions of (Q), the random term rf is represented by 
choosing a random number in every iteration, uniformly 
distributed in the range — 0.5ry < rf < O.brj. In the follow- 
ing we use the same choice of parameters as in [ fl6[ with 
A = 3m/ sec 2 , and a numerical time interval At = O.lsec. 

(a) Signalized intersection: In an oversaturated 
signalized intersection, the congested flow in each direc- 
tion is affected only by the parameters of the traffic sig- 
nals of this direction. Therefore we consider a single di- 
rection for simplicity (as e.g. in |L9||). To study this 
case, a single traffic light is placed at position L/2 on 
a road with length L and periodic boundary conditions. 
The flux / is measured at the position of the traffic light. 
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Simulations are performed with different durations of the 
red, amber (yellow), and green lights (r r , r y [25f[ , r g re- 
spectively), and of r_ which is the total time in each 
signal period where the intersection is not in use by any 
of its approaches. Therefore r g = t — t_ — (T r + T y ) where 
r is the total signal period. P r = t t /t is the relative du- 
ration of the red light. 

Our goal here is to optimize the flux in a specific road 
without affecting the parameters of the whole complex 
system of roads. We therefore consider r_, t v and espe- 
cially P r as given constraints, and explore the relation 
between the flux / and the signal period t. Typical ex- 
amples of this relation are shown in Fig. 1. 

Apart from the trivial flux oscillations explained below, 
Fig. 1 shows a monotonic increase of the flux as r grows, 
for the deterministic model (?/ = 0). For the stochastic 
model (i] > 0), however, an optimal signal period can be 
clearly seen. In this case, a crossover is observed with in- 
creasing r from the deterministic /(r) to decreased values 
of flux. The optimal value of r approaches the crossover 
point for small rj. 



o 

CD 



LL 



0.24 - 



0.22 




0.18 
0.25 



1 2 4 
Signal period i [min] 



16 



/ = fa 



foil- Pr 



T- 
T 



(2) 



where /o does not depend on r (see jjlj and references 
within) . For the signal period values r t_ that are dis- 
played in Fig. 1, one would expect to find /(r) « const. 



according to (g). Therefore the numerical finding of the 
influence of r on the flux shown in Fig. 1, revealing an 
optimal signal period, requires deeper investigation. 

In an attempt to explain our findings, we performed 
extensive simulations of the deterministic model Eq. (|l|) 
on homogeneous systems with periodic boundary con- 
ditions, starting from different initial conditions |p6| . 
Fig. 2a presents the flux measured for stable stop-and-go 
(surface) and homogeneous (thick lines) states, showing 
our finding that the traditional 'fundamental diagram' 
(density-flux relation) has to be generalized into a three- 
dimensional density-wavelength-flux relation. 

The projection of the surface in Fig. 2a on the density- 
flux plane (thin curves in Fig. 2a) provides a two- 
dimensional region in this plane, qualitatively similar to 
that found empirically by Kerner for synchronized 
flow. In our system this phenomena occurs since the 
flux of the periodic states depends not only on the den- 
sity but also on the state's wavelength A p7| . That is, 
for a given value of density, there exist many stable pe- 
riodic states which have a range of wavelengths in the 
density-flux plane. Since the flux is also determined by 
the wavelength (see Fig 2b), the density-flux relation be- 
comes multi- valued. 




FIG. 1. Relation between signal period and flux for the val- 
ues of acceleration noise amplitude r\ — 0, 0.5, 2, 5, 10m/ sec 2 

(top to bottom). Traffic lights parameters are P r = 1/3, a > 

t v — t_ = 2sec. The total number of cars in the system is " 3 

N — 400 and its length is L = lOfcm. The nine open circles § 
correspond to the nine instances presented in Fig. 3. 

The common traffic flow and signal control theories do 
not predict such influence of r on the flux itself (see e.g. 
| [18| ^0|,|3| ) . The average flow in each direction of a sig- 
nalized intersection /(r) for our consideration is expected 
to be 
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FIG. 2. (a) Density-wavelength-flux relation for the differ- 
ent states of the deterministic model, including stable peri- 
odic states (surface), stable and unstable homogeneous states 
(thick and dotted lines, respectively). Some curves with fixed 
wavelength (A" 1 = 2/60, 3/60, 12/60) were projected on 
the density-flux plane (thin curves), (b) A cross-section of 

(a) for a density p = 0.06veh/m, demonstrating the typical 
dependence of the flux on the wavelength, (c) Noise stability 
threshold amplitude rjth, above which the states presented in 

(b) become unstable. 

Fig. 2c shows the noise stability threshold amplitude 
rjth, above which the states presented in Fig. 2b become 
unstable. A comparison of the last two figures shows 
that states with flux values close to the minimal are the 
most stable in the presence of noise, while states with 
higher values of flux (close to the upper bound of the 
two dimensional region) are metastablc. Therefore real- 
life stop-and-go traffic might be expected to show values 
of flux far below the optimum. 

This new insight on the nature of the two dimensional 
region in the density-flux plane provides also an explana- 
tion to our finding of a non-trivial relation between flux 
and signal period. For rj = 0, the monotonically increas- 
ing /(t) displayed in Fig. 1 corresponds to the left branch 
of Fig. 2b. High values of r stimulate transition to states 
with high values of A - and therefore with high values of 
flux, according to Fig. 2b. These metastable states sur- 
vive since there is no noise ^8|. The oscillations in flux 
observed for relatively small r's can be easily explained 
& 

The crossover observed in Fig. 1 for the stochastic 
model (ry > curves) can be related to a crossing of the 
noise stability threshold rjth- For values of r below the 
crossover point, the value of /(t) is close to that of the de- 
terministic model, since rj < r\ t h- When r (and therefore 
A) is increased, two trends are expected according to Figs 
2b and 2c: an increase in the flux and a decrease in rjth-, 
until 77 > r/th, where a crossover to values of flux lower 
than that of the deterministic case. Therefore the opti- 
mal t is usually close to the crossover point. Note that 
this crossover from deterministic to non-deterministic be- 
havior with changing t occurs in spite of the fact that the 
noise amplitude is fixed, and is related to the r (and A) 
dependence of rjth- 

We therefore see that the non-trivial flux-wavelength 
relation is the reason for the unexpected behavior of /(t). 
The deviations between the theoretical prediction of (g) 
and the numerical measurements presented in Fig. 1 and 
the crossover observed in f(r) for 77 > are related to 
the relation between flux and wavelength discussed above 
and to differences in noise stability threshold of different 
periodic states. 

To visualize the effect of the signal period and the ac- 
celeration noise amplitude on the flow, nine space-time 
diagrams are presented in Fig. 3a, below, at, and above 
the crossover. These nine diagrams correspond to the 
nine instances denoted by open circles in Fig. 1. Each 



dot in these space-time diagrams represents the position 
of a single car at a certain time. The dark regions there- 
fore show the dense regions on the road. Looking at the 
diagrams in Fig. 3a, one can see the increase of the dom- 
inant wavelength with increasing r. But unlike the de- 
terministic case where the flow is periodic and the flux is 
increasing with A and therefore with r, in the stochastic 
model (rj > 0) small jams emerge in the low density re- 
gions when the noise or the signal period exceed certain 
thresholds. In such cases, other values of wavelength, 
smaller than that induced by the traffic light are effec- 
tively involved, resulting values of flux lower than that of 
the X] = case (see Fig. 1). 

For further support of this interpretation we evaluate 
the periodicity in the flow using single vehicle data col- 
lected at the intersection. Inspired by we calculate 
the auto-covariance ac v (t) of the velocity function v(t') 
measured at the intersection, 



ac v (t) 



< v(t')v(t' + t)> - < v(t') >< v(t' + t)> 



< v(t') 2 > - < v{t') > 2 



(3) 



where a linear continuation of the discrete function v(t') 
is used. The brackets < . . . > indicate averaging over 
time t' . Displayed in Fig. 3b are the auto-covariance 
functions for the 9 instances of Fig. 3a, respectively. As 
can be seen from this figure, ac(t — t) — 1 for all the 
i] = instances, implying the flow for these cases is com- 
pletely periodic, and that the time period is equal to the 
period of the traffic light. For r\ > 0, ac v (t — r) decreases 
as rj or r are increased. This decrease is related to the 
appearance of small jams as can be seen in Fig. 3a. The 
decrease in ac v (t — t) with r also implies that noise be- 
comes more effective as t grows. A comparison of Fig. 
3b to Fig. 1 shows that when ac v (t = r) « 1, the flux 
approaches the deterministic value. 
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FIG. 3. (a) Space-time diagrams and (b) auto-covariance 
functions of systems with single traffic light with parameters 
as the nine instances denoted with open circles in Fig. 1. The 
position of the traffic light is at x — 5km. Gray lines in (b) 
corresponds to 77 = 0, solid lines to rj = 2, and dashed lines 
to j] = 5. 
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FIG. 4. (a) Relation between on-ramp signal period and 
flux on the main road for P g — 0.3,0.5,0.6,0.8,1.0 (top to 
bottom) and noise amplitude rj = 2m/ s 2 . The total number 
of cars in the system is N — 300, its length is L = 10km, and 
the flux is locally measured on the main road 100m upstream 
to the on-ramp. Here f in = 0.1 and f ma x = 0.333. (b) A 
comparison between the optimal flux (upper curve) and the 
flux without the presence of a traffic light (lower curve), as a 
function of r\. 



(b) On-ramp with signalized entrance: Next, 
we study another example of a localized inhomogeneity, 
caused by an on-ramp. To make this inhomogeneity peri- 
odic, we introduce traffic signals at the downstream end 
of the ramp, and study its effect on the flux in the main 
road. We focus on cases where the average incoming flux 
fin is high enough to induce congestion on the main road 
(see H), but still low enough to avoid congestion on the 
secondary road. Similarly to |l7|] , we introduce also an 
off-ramp with equal flux of outgoing vehicles f ou t = f% n 
in a large distance from the on-ramp, so that the total 
number of cars in the system is conserved. The entrance 
and the exit of cars from the ramps are performed in an 
adiabatic manner as in |l6[ |. The incoming vehicles are 
allowed to enter the main road during the green light 
period r g , and are stopped during the red light period 
T r . Here the signal period is r = r g + r r and the rela- 
tive duration of the green light is P g = t 9 /t. But unlike 
the signalized intersection where P r was predetermined, 
here P g is one of the optimization parameters, in addi- 
tion to t. The range of possible values for this parameter 
is fin/ fmax < P g < 1 where f max is the maximal possi- 
ble value of the incoming flux. This lower bound of P g 
is considered to avoid congestion on the secondary road, 
since cars approach the queue near the traffic light with 
rate fi n , and this queue is discharged with rate f max dur- 
ing green light. The upper bound P g — 1 is related to an 
unsignalized on-ramp. 



Typical relations between flux and signal period are 
shown in Fig. 4a, for different values of P g and for 
i] = 2m/s 2 . The lower curve corresponds to an unsignal- 
ized on-ramp (P g = 1), so it is found that the intro- 
duction of a traffic light increases the flux on the main 
road. Since the lower bound of P g ensures that the av- 
erage influx from the on-ramp remain the same as in the 
unsignalized on-ramp, we can say that the increases in 
the flux on the main road is obtained without causing 
congestion on the secondary road. 

As can be seen from Fig. 4a, for each value of P g there 
usually exist a single maximum in the flux. The short 
relative durations of green light P g « fi n / f max usually 
(but not always) yields a higher flux. A comparison of 
the maximal flux at the optimal signal parameters to the 
flux without a traffic light is presented in Fig 4b, for dif- 
ferent values of noise. The relative increase in the flux 
due to the introduction of a traffic light varies from 1.0% 
for rj = 10m/s 2 , through 10.0% for rj = 2m/s 2 up to 
13.9% for rj = 0. 
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